plasma Module

модуль параметров плазмы



Variables

Type Visibility Attributes Name Initial
integer, public :: ngrid

ASTRA radial grid number

integer, public :: nspl

ASTRA radial grid number

real(kind=wp), public :: tcur

время (придумать название для переменной получше)

real(kind=wp), public :: rm

minor radius in mid-plane, cm

real(kind=wp), public :: b_tor0

тороидальное магнитное поле временно нужно две переменных, тоже нужно исправить

real(kind=wp), public :: b_tor

тороидальное магнитное поле

real(kind=wp), public :: r0
real(kind=wp), public :: z0
real(kind=wp), public :: rh1
real(kind=wp), public, dimension(:), allocatable :: con
real(kind=wp), public, dimension(:), allocatable :: tem
real(kind=wp), public, dimension(:), allocatable :: temi
real(kind=wp), public, dimension(:), allocatable :: zeff
real(kind=wp), public, dimension(:), allocatable :: afld
real(kind=wp), public, dimension(:), allocatable :: rh
real(kind=wp), public, dimension(:), allocatable :: rha
real(kind=wp), public, dimension(:), allocatable :: drhodr
real(kind=wp), public, dimension(:), allocatable :: delta
real(kind=wp), public, dimension(:), allocatable :: ell
real(kind=wp), public, dimension(:), allocatable :: gamm
real(kind=wp), public, dimension(:), allocatable :: amy
real(kind=wp), public :: tet1

common /a0a2/

real(kind=wp), public :: tet2

common /a0a2/

real(kind=wp), public :: xmi

/a0ef3/ xmi,cnye,cnyi,xsz,vt0

real(kind=wp), public :: cnye

/a0ef3/ xmi,cnye,cnyi,xsz,vt0

real(kind=wp), public :: cnyi

/a0ef3/ xmi,cnye,cnyi,xsz,vt0

real(kind=wp), public :: xsz

/a0ef3/ xmi,cnye,cnyi,xsz,vt0

real(kind=wp), public :: vt0

/a0ef3/ xmi,cnye,cnyi,xsz,vt0

real(kind=wp), public :: cnstvc
real(kind=wp), public :: ww

common /a0ef2/ ww

real(kind=wp), public :: cltn

common /a0ef1/ cltn

real(kind=wp), public :: vperp(50,100)
real(kind=wp), public :: cnstal
real(kind=wp), public :: zza
real(kind=wp), public :: zze
real(kind=wp), public :: valfa
real(kind=wp), public :: vpmax
real(kind=wp), public :: vk(100)
real(kind=wp), public :: sk(100)
integer, public, parameter :: ipsy = 5

ipsy = number of polinomial decomposition coefficients used for interpolation of Zakharov's moments.

integer, public, parameter :: ncoef = 5

ipsy = number of polinomial decomposition coefficients used for interpolation of Zakharov's moments.

real(kind=wp), public, dimension(ipsy) :: cdl
real(kind=wp), public, dimension(ipsy) :: cly
real(kind=wp), public, dimension(ipsy) :: cgm
real(kind=wp), public, dimension(ipsy) :: cmy
real(kind=wp), public, dimension(ipsy) :: coeffs
real(kind=wp), public :: y2dn(501)

common /a0l3/

real(kind=wp), public :: y2tm(501)

common /a0l3/

real(kind=wp), public :: y2tmi(501)

common /a0l3/

real(kind=wp), public :: y2zeff(501)

common /a0l5/

integer, public :: ncheb
real(kind=wp), public :: chebne(50)

common/ne_cheb

real(kind=wp), public :: chebdne(50)

common/ne_cheb

real(kind=wp), public :: chebddne(50)

common/ne_cheb

real(kind=wp), public :: enorm(100)

em поле и еще что-то

real(kind=wp), public :: fst(100)

em поле и еще что-то


Functions

public function fn(x)

plasma density, cm^-3

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: x

Return Value real(kind=wp)

public function fvt(r)

нет описания

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: r

Return Value real(kind=wp)

public function fn1(x, fnp)

plasma density and its derivative

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: x
real(kind=wp), intent(out) :: fnp

Return Value real(kind=wp)

public function fn2(r, fnp, fnpp)

plasma density and its first and second derivatives

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: r
real(kind=wp), intent(out) :: fnp
real(kind=wp), intent(out) :: fnpp

Return Value real(kind=wp)

public function ft(x)

electron temperature, erg ft=y ! kev

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: x

Return Value real(kind=wp)

public function fti(x)

ion temperature, kev

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: x

Return Value real(kind=wp)

public function zefff(x)

z_effective profile

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: x

Return Value real(kind=wp)

public function obeom(ptet, pa)

Arguments

Type IntentOptional Attributes Name
real :: ptet
real :: pa

Return Value real(kind=wp)

public function ploshad(ptet, pa)

Arguments

Type IntentOptional Attributes Name
real :: ptet
real :: pa

Return Value real(kind=wp)

public function gaussint(f, a, b, r, eps)

что-то про гаусс

Arguments

Type IntentOptional Attributes Name
real(kind=wp) :: f
real(kind=wp) :: a
real(kind=wp) :: b
real(kind=wp) :: r
real(kind=wp) :: eps

Return Value real(kind=wp)


Subroutines

public subroutine init_plasma(NA1, ABC, BTOR, RTOR, UPDWN, GP2, AMETR, RHO, SHIF, ELON, TRIA, MU, NE, TE, TI, ZEF, UPL)

variant amy=(btor/q)rho(drho/dr) is a function of "minor radius" r=rh(i). Poloidal magnetic field: B_pol=amy(r)*sqrt(g22/g), where g is determinant of 3D metric tensor and g22 is the (22) element of the tensor, normalized on ABC^4 and ABC^2, correspondingly.

Read more…

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: NA1
real(kind=wp), intent(in) :: ABC
real(kind=wp), intent(in) :: BTOR
real(kind=wp), intent(in) :: RTOR
real(kind=wp), intent(in) :: UPDWN
real(kind=wp), intent(in) :: GP2
real(kind=wp), dimension(*) :: AMETR
real(kind=wp), dimension(*) :: RHO
real(kind=wp), dimension(*) :: SHIF
real(kind=wp), dimension(*) :: ELON
real(kind=wp), dimension(*) :: TRIA
real(kind=wp), dimension(*) :: MU
real(kind=wp), dimension(*) :: NE
real(kind=wp), dimension(*) :: TE
real(kind=wp), dimension(*) :: TI
real(kind=wp), dimension(*) :: ZEF
real(kind=wp), dimension(*) :: UPL

public subroutine init_parameters()

!!!!!! ptkev=ft(zero)/0.16d-8 !Te in keV ccurnr=pqepqe0.333d-9/pme

Arguments

None

public subroutine write_plasma(time_stamp)

Arguments

Type IntentOptional Attributes Name
real(kind=wp), intent(in) :: time_stamp

public subroutine write_lcms()

write lcms

Arguments

None

public subroutine find_volums_and_surfaces()

Arguments

None

public subroutine calc_enorm()

fst(j)=pnxlogjc04/pi4/vt3

Arguments

None

public subroutine init_maxwell()

Arguments

None